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ABSTRACT 

Recent observations of high energy cosmic-ray electrons by the Fermi-LAT and the HESS 
experiments between 20 GeV and 5 TeV have found that the energy spectrum closely follow a 
broken power-law with a break at around 1 TeV. On the other hand, measurements of cosmic- 
ray secondary-to-primary ratios like the boron-to-carbon ratio seem to indicate a possible 
change in the slope at energies around 100 GeV/n. In this paper, we discuss one possible 
explanation for the observed break in the electron spectrum and its possible correlation with 
the flattening in the secondary-to-primary ratios at higher energies. In our model, we assume 
that cosmic -rays after acceleration by supernova remnant shock waves, escape downstream 
of the shock and remain confined within the remnant until the shock slows down. During 
this time, the high-energy electrons suffer from radiative energy losses and the cosmic-ray 
nuclei undergo nuclear fragmentations due to their interactions with the matter. Once the 
cosmic -rays are released from the supernova remnants, they follow diffusive propagation in 
the Galaxy where they further suffer from radiative or fragmentation losses. 
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1 INTRODUCTION 

Few years ago, measurements of high energy cosmic-ray (CR) elec- 
trons above 10 GeV were performed mainly by balloon-borne ex- 
periments (Kobayashi et al. 2004 and references therein). Their 
measurements showed that the electron spectrum follow a power- 
law behaviour of the form E~ r with the index T « 3.2 without any 
significant features up to energies around 2 TeV. However, recent 
measurements made by the ATIC balloon experiment had found a 
sharp peak at E « 600GeV (Chang et al. 2008). But, this feature 
was not later detected by two other experiments, the spaced-bome 
Fermi-LAT and the ground-based HESS experiments (Abdo et al. 
2009, Aharonian et al. 2008b, 2009) which measured the energy 
spectrum in the range of 20GeV — 5TeV. Their combined spec- 
trum can be closely represented by a broken power-law with spec- 
tral indices Ti ~ 3 and T2 ~ 4 for energies below and above 1 
TeV respectively. In standard CR propagation studies assuming a 
homogeneous source distribution, one way to explain such a spec- 
tral behavior is to assume that the electron source spectrum fol- 
low a power-law behavior with an exponential cut-off of the form 
exp(— E/E c ) with E c m ITeV. Radio and X-ray observations do 
support such a form of electron spectrum inside supernova rem- 
nants (SNRs), but with cut-off energies as high as E c ~ 80 TeV as 
suggested by the study of Reynolds & Keohane 1999 considering 
magnetic field strengths of 10/iG inside the remnants. However, for 
young SNRs where magnetic field amplification seems to occur, the 
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field strengths can reach values even more than 100/iG (Volk et al. 
2005) and the maximum energies for the electrons can be strongly 
limited by radiative losses to values less than ~ 10 TeV. 

The observed break can also be an effect of inhomogeneous 
distribution of CR sources contributing to the electrons at higher 
energies (Atoyan et al. 1995, Kobayashi et al. 2004, Blasi 2009). 
This is because high energy electrons cannot travel far distances 
in the Galaxy due to their faster radiative cooling. Therefore, it is 
possible that most of the TeV electrons that we measure in the So- 
lar System are produced by few young nearby sources. Then, the 
general assumption of a continuous source distribution may break 
down and a more reasonable treatment would be to take into ac- 
count a discrete source distribution. A common method to do this 
is to use the idea of separation of distant and nearby sources, i.e., 
assuming a continuous distribution for the distant sources and con- 
sidering known SNRs or pulsars as the nearby discrete sources. One 
of the main uncertainties involved in such studies can be the effect 
of missing sources due to detection biases. Moreover, even for the 
known sources, the lack of precise informations about the source 
parameters like the distance, age and the spectral informations can 
lead to strong uncertainties in the high energy electron spectrum. 
Recently Delahaye et al. 2010 presented a detailed study using this 
method where they included all the known SNRs and pulsars lo- 
cated within 2 kpc from the Earth and they found that the total high 
energy spectrum from the nearby sources strongly depends on sev- 
eral source parameters including the assumed cut-off energy. They 
also showed that the cut-off energy should be somewhere around 
a few TeV in order to explain the break in the observed data. An 
alternative explanation for the break, as will be shown in this pa- 
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per, is that after acceleration the high energy electrons might have 
suffered significant losses within the sources itself before they are 
released into the Galaxy. 

On the other hand, measurements of CR secondary-to-primary 
(s/p) ratios like the boron-to-carbon (B/C) up to E « 100 GeV/n 
by several independent experiments (see the experiments listed in 
Stephens & Streitmatter 1998) show that above 1 GeV/n, the ratio 
decreases with energy as ~ E~ 5 with the index 8 « 0.6 imply- 
ing an energy dependent CR propagation path length in the Galaxy. 
However, these measurements seem to indicate a possible harden- 
ing at higher energies (see also Swordy et al. 1990). Recent data 
from the CREAM balloon-borne experiment also show flattening 
at energies around 1 TeV/n (Ahn et al. 2008). One way to explain 
this is to assume that CRs traverse some minimum amount of mat- 
ter of approximately 0.3 g cm~ 2 in the Galaxy (Ave et al. 2009). 
This implies an energy dependent escape path length of CRs up 
to some certain energy and beyond that energy, the path length 
becomes constant in energy. This view is, in fact, supported by 
the observed CR anisotropy which remain almost constant above 
around 100 GeV energies. But such model need to assume a break 
in the source spectrum in order to explain the observed CR spec- 
trum which closely follow a pure power-law without any break up 
to the knee. Another possible explanation is that some amount of 
secondaries might be produced inside the sources due to the nuclear 
interaction of the primaries with the matter (Berezhko et al. 2003). 
Under such models, the CR anisotropy is expected to increase with 
energy as E s which does not agree well with the observed data. 
However, it is quite possible that at higher energies the anisotropy 
might not be determined by the global diffusion leakage of CRs 
from the Galaxy. Rather, it might be the effects of high energy CRs 
coming from nearby sources because of their faster diffusion in 
the Galaxy (Ptuskin et al. 2005, Thoudam 2007). One more pos- 
sible explanation which we will not discuss here is the possible 
reacceleration of some fraction of the background CRs by strong 
SNR shocks while propagating in the Galaxy (Wandel et al. 1987, 
Berezhko et al. 2003). 

In this paper, we present one possible explanation for the ob- 
served break in the electron spectrum at around 1 TeV and its possi- 
ble correlation with the flattening in the B/C ratio at energies above 
around 100 GeV/n. In our model, we assume that CRs are accel- 
erated by SNR shock waves by diffusive shock acceleration (DSA) 
mechanism (Bell 1978, Blandford & Eichler 1987). Under DSA 
theory, suprathermal particles in the tail of the Maxwell-Boltzmann 
distribution present in the interstellar medium (ISM) are injected 
into the supernova shock front. They are then reflected back and 
forth several times across the shock by the magnetic turbulence 
generated on either side of the shock and in each crossing, the par- 
ticles gain energy by first order Fermi acceleration. In a simple pla- 
nar shocks model, such a mechanism naturally leads to a power-law 
spectrum of the form E~ F with F = 2 for strong shocks which is 
in good agreement with radio observations of several SNRs (Green 
2009). 

One important consideration of the standard DSA theory is 
that most of the particles do not escape upstream of the shock, they 
can only escape downstream mainly due to advection by the bulk 
flow and remain confined within the remnant due to the strong mag- 
netic turbulence generated by the CRs themselves. Though it is still 
not fully understood, it is generally considered that the particles 
remain inside the remnant until the shock remains strong. During 
the confinement period, the high energy electrons may suffer from 
radiative energy losses while the nuclear components suffer from 
nuclear fragmentations. At later stages when the shock slows down 



and does not efficiently accelerate the CRs, the magnetic turbulence 
level goes down and the particles can no longer be confined effec- 
tively within the remnant. At such stage, all the particles present 
inside the remnant escape and they are injected into the ISM. For a 
typical ISM density of 1 H cm - ' , this happens at around 10 yr af- 
ter the supernova explosion (Berezhko & Volk 2000). The scenario 
described above might be true mostly for the lower energy particles 
because it is quite possible that some of the highest energy particles 
may escape the remnant already at the start of the sedov phase it- 
self due to their faster diffusion. Such an energy dependent escape 
of particles has been discussed in a number of literatures (Ptuskin 
& Zirakashvili 2005, Gabici, Aharonian & Casanova 2009). Recent 
high energy 7-ray observations of SNRs associated with molecular 
clouds also suggest that some of the high energy particles might 
have already escaped the remnant at much early times (Aharonian 
et al. 2007, Aharonian et al. 2008a). But, as mentioned above, it is 
still not clear how and when exactly the particles escape the rem- 
nant. For simplicity, we consider an energy independent scenario 
for our present work and assume that all the particles are released 
into the ISM at some characteristic time after the supernova explo- 
sion. 

Once the CRs are released into the ISM, we assume that they 
undergo diffusive propagation in the Galaxy where the electrons 
again suffer from radiative losses and the nuclei from nuclear spal- 
lation with the interstellar matter. Therefore, in our model the CR 
spectra that we finally observe in the Solar System are modified 
from their original source spectrum (the one generated by the SNR 
shocks) due to the various interactions or energy loss processes oc- 
curring not only during their propagation in the Galaxy, but also 
within the SNRs itself. 

Our paper is planned as follows. In section 2, we present our 
calculations for the CR spectra inside the SNRs and in section 3, we 
calculate the spectra in the Galaxy. Then, in section 4, we compare 
our results with the observed data and also give a short discussions 
about our results. Finally in section 5, we give a brief conclusion of 
our study. 



2 CR SPECTRA WITHIN THE SNRS 

We assume that CR acceleration by SNR shock waves begins at 
the time of the supernova explosion itself and the CRs then escape 
downstream of the shock and remain confined within the remnant. 
We further assume that at some later stage of the SNR evolution 
characterize by age t = T when the shock slows down, particle 
acceleration stops and also all the CRs present inside the remnant 
are released into the ISM. In our study, we do not consider the 
expansion of the remnant and hence neglect the effect of adiabatic 
energy losses and other related effects. Detailed study including 
the evolution of the remnant, the weakening of the shock with time 
and the possible energy dependent escape of CRs will be presented 
elsewhere. 



2.1 High energy electrons 

Under the assumption that the acceleration time of CRs is much 
less than their confinement time within the SNR, the CR electron 
spectrum inside an SNR can be described by the following equa- 
tion, 



(1) 



A possible correlation between the high energy electron spectrum and the cosmic-ray secondary-to-primary ratios 



where N e (E, t) is the total number of electrons of kinetic energy 
E present within the SNR at time t, b(E) — —dE/dt is the energy 
loss rate and q e (E) is the source term which, in our case, is the 
rate at which CR electrons are injected downstream of the shock. 
We believe that Eq. (1) represents a valid approximation at least 
for energies up to around 10 TeV for which the typical acceleration 
timescale is less than ~ 100 yr (Sturner et al. 1997). In Eq. (1), we 
consider the energy loss of the electrons to be due to synchrotron 
and inverse compton interactions which is true for energies E > 10 
GeV. Therefore, we take in the Thompson regime 



where, 



b(E) = aE 2 



(2) 



with a = 1.01 x 10~ 16 (w p h + wb) GeV -1 s _1 , where w p h and 
wb are the energy densities of the radiation fields and the magnetic 
field respectively in eV cm -3 . With this, we define the radiative en- 
ergy loss timescale for the electrons as ti oas (E) = l/(aE). Now, 
for the source spectrum taken as power-law of the form, 



q e (E) = k e E~ 



(3) 



the solution of Eq. (1) at time t can be written as (Kardashev 1962, 
Gratton 1972) 

ME,t) = ^ [l - (1 - aEtf- 1 ^ , for E < E 



(T-l)o 
fc e g-< r+1 > 
(T-l)o 



for E^E (4) 



where Eq = l/(at) is the energy at which the energy loss time 
ti OS3 becomes equal to time t. Eq. (4) shows that for energies E -C 
Eq, the electron spectrum still reflects the source spectrum because 
of their large ti oss , i.e. N e (E) oc E~ F . However, for electrons with 
energies £» Eq, their spectrum is quite steeper due to their faster 
energy loss rate and follows N e (E) oc i5~( r+1) . 



2.2 CR nuclei 

For the CR primary nuclei, the spectrum inside the SNR can be 
described by 



dN, 
dt 



- + ij'capNp = q p 



(5) 



where N P (E, t) is the total number of particles of kinetic energy 
per nucleon E at time t, rj is the matter density in the SNR, c is 
the velocity of light and a p is the spallation cross-section of the 
primary nuclei assumed to be independent of energy. In Eq. (5), we 
assume that the source term q P (E) = k p E~ F has the same index 
as those of the electrons. 

The solution of Eq. (5) can be obtained as 



il'cGp L 



e -i7' Ct Tptj 



(6) 



For time much less than the nuclear spallation time, i.e. t <g; 
l/(r/ cop), Eq. (6) becomes N P (E, t) « q p (E)t. 

During the time when primary CRs are confined within the 
SNRs, they interact with matter and produce secondary nuclei of al- 
most the same kinetic energy per nucleon as their primaries. These 
secondaries can also be described by an equation similar to Eq. (5) 
by replacing the source term by rj ca ps N p (E, t) where a ps is the 
total fragmentation crossection of primary to secondary. The solu- 
tion for these secondaries is then obtained by 



NJE,t) = [i_ e -iW + F ] 

cr p r)'cG s L J 



(7) 



(cr s 



In Eq. (7), a s represents the spallation cross-section of the sec- 
ondary nuclei. Now, taking t = T which is the CR confinement 
time inside the SNR, we can use Eqs. (4), (6) & (7) to calculate 
the spectrum of CR electrons and the spectra of primary and sec- 
ondary CR nuclei finally injected into the ISM from a single SNR. 
Then, knowing the rate of supernova explosion per unit volume in 
the Galaxy, we can calculate the rate at which CRs are injected per 
unit volume in the Galaxy. 



3 CR SPECTRA IN THE GALAXY 

CRs after escaping from the SNRs undergo diffusive propagation 
in the Galaxy due to scattering either by magnetic field irregu- 
larities or by self excited Alfven and hydromagnetic waves. For 
the present work, we assume the diffusion region as a cylindri- 
cal disk of infinite radius with finite half-thickness H and that the 
sources as well as the matter are distributed uniformly and con- 
tinuously in the Galactic disk with half-thickness h and radius R, 
where R ^> H ^> h. The details of the geometry are described in 
Thoudam 2008. During the propagation, high energy electrons in- 
teract with the background radiation and the magnetic fields and 
lose their energies. On the other hand, CR nuclei undergo nu- 
clear spallation interactions with the interstellar matter and produce 
lighter nuclear species. 

As already mentioned in section 1, the assumption of a con- 
tinuous source distribution may not be fully appropriate for high 
energy electrons particularly those in the TeV region because of 
their faster energy loss rate. Electrons with energies greater than 
1 TeV cannot travel distances more than ~ 1 kpc in the Galaxy 
through diffusive propagation before they lost all their energies. 
Therefore, high energy electrons from distant and old sources may 
not reach the Earth effectively and TeV electrons that we observe 
can be mostly dominated by those produced by few young local 
sources. The effect of this can be that the spectrum at high ener- 
gies can be quite complex because of its strong dependence on the 
local source parameters (Delahaye et al. 2010). Moreover, it may 
even show up features related to the stochastic nature of the local 
sources in space and time (see e.g., Pohl & Esposito 1998). 

In our study, our main focus is to explore the possibility of ex- 
plaining the break in the electron spectrum as an effect of CR con- 
finement in the downstream region of SNRs and at the same time, 
following the standard model of CR propagation in the Galaxy. 
Such an effect of confinement within the sources are generally not 
considered in CR propagation studies and we believe that they may 
exist if SNRs are the main sources of galactic CRs. Therefore, we 
do not intend to focus on a detailed source distribution and in what 
follows, we adopt the continuous and stationary source distribution 
for calculating the CR spectrum in the Galaxy. 



3.1 High energy electrons 

Under the diffusion model, the propagation of high energy elec- 
trons in the Galaxy can be described by 



V • (DVn e ) + JL{b{E)n e } = -Q e 



(8) 



where n e (r, E) is the electron number density in the Galaxy, 
D(E) is the diffusion coefficient which is assumed to be constant 
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through out the Galaxy and Q e (r, E) is the source term given by 
Q e (r, E) = $tN e (E, T)S(z). Here, N e (E, T) is given by Eq. (4) 
which represents the total amount of CR electrons liberated by an 
SNR and Ji denotes the rate of supernova explosion per unit surface 
area on the Galactic disk. As already mentioned, high energy elec- 
trons cannot travel large distances because of their large energy loss 
rate. Therefore, their equilibrium spectrum is not much affected by 
the presence of the halo boundary H and we solve Eq. (8) without 
imposing the finite boundary conditions for our study. The solution 
at the position of the Earth (z = 0) is then obtained as 

-i -1/2 

D(u) 



n e (E) = 



2^b(E) 



f 

J E 



dE'N e (E',T) 



I 

J E 



b{u) 



du 



(9) 



For the diffusion coefficient D(E) oc E s , Eq. (9) shows that 



n e (E) oc 



b{E) 



dE'- 



N e (E',T) 
yjE 6 - 1 - E' s - 



(10) 



From Eq. (10), we see that for a power-law source spectrum such 



N e (E) 



oc E 7c , we get n e (E) oc E 



-1+/3 



where /3 

1 — 5)/2. Thus, the equilibrium electron spectrum in the Galaxy 
is modified mainly by the radiative energy losses with a small cor- 
rection P arising due to the diffusive propagation effect. 



3.2 CR nuclei 

The steady state transport equation for primary CR nuclear species 
in the Galaxy, neglecting the effect of energy losses due to ioniza- 
tion, is given by 

V • (DX7n p ) - 2hr)ca p 8(z)n p = -Q p (11) 

where n p (r, E) is the density of primary CRs of energy per nu- 
cleon E, r\ is the matter density in the ISM and Q p (r,E) — 
ftN p (E,T)S(z) is the source term where N P (E,T) is given by 
Eq. (6). Eq. (11) is solved by imposing proper boundary conditions 
as described in detail in Thoudam 2008 and the solution at z = 
is given by 

FvRN p (E, T) f°° J 1 (KR)dK 



n p (E) 



2D 



f 

Jo 



[Kcoth(KH) + '^p-] 



(12) 



where Ji is the Bessel function of order 1. 

Primary nuclei during their propagation in the Galaxy interact 
with matter in the Galaxy and produce secondary nuclei. The equi- 
librium spectrum for these secondaries n s (r, E) in the Galaxy can 
be described by 



V • (DVn s ) — 2hrjca s 5(z)n s 



-Qs 



(13) 



where the source term Q s (r, E) — 2hr\ca ps n p (r, E)S(z). There is 
also an additional component of the secondaries in the Galaxy due 
to their production inside the SNRs as discussed in section 2.2. The 
equilibrium spectrum n' a for this additional component can also be 
described by an equation similar to Eq. (13) with the source re- 
placing by Q' s (r, E) = UN 3 (E, T)S(z), where N s (E, T) is given 
by Eq. (7). The spectra of these two secondary components at the 
position of the Earth can be obtained respectively as 

( v\ o*, tm R f°° Ji(KR)dK 

n s (E) = 2hr)ccj 1 ~ ' '• 



[Kcoth(KH) + 2^ 



and 



n's(E) = 



RUN a (E,T) 
2D 



f 

Jo 



J 1 {KR)dK 



[Kcoth(KH) + 



2hrjc 



(14) 



(15) 



Having calculated the spectra for the primaries and the secon- 
daries in the Galaxy, the overall s/p ratio can be obtained by simply 
taking the ratio (n 3 + n' s )/n p . Note that in the standard model of 
CR propagation in the Galaxy where secondary nuclei are assumed 
to be produced only in the ISM, the ratio is given by n s /n p . 



4 RESULTS AND DISCUSSIONS 

For our calculations, we choose R = 16 kpc, H = 5 kpc, h — 200 
pc and the rate of supernova explosion as Jft = 25 Myr -1 kpc -2 
(Grenier 2000). This corresponds to a total supernova explosion 
rate of ~ 1/50 yr in our Galaxy. The total energy density of the 
radiation fields is taken as w ph — wmbr + w op + wfir where 
wmbr, w op and wfir are the energy densities of the microwave 
background, the NIR-optical and the FIR radiation fields respec- 
tively. We take wmbr = 0.25 eV cm~ 3 , w op — 0.5 eV cm -3 
(Mathis, Mezger & Panagia 1983), w FIR = 0.2 eV cirT 3 (Chi & 
Wolfendale 1991) and the magnetic field value as B — 6/iG (Beck 
2001). We assume these values to be same for both in the Galaxy 
and in the SNRs. We take the ISM density as r\ = 1 cm -3 (Gordon 
& Burton 1976, Bronfman 1988) typical for our Galaxy and the CR 
diffusion coefficient in the Galaxy as (Thoudam 2008) 

D(E) = D (^) , for£<^ 



A, ( | ) - to, L >L 



(16) 



where D = 2.9 x 10 28 cm 2 s _1 and S = 0.6. For CR nuclear 
species, E represents energy per nucleon and the particle rigidity 
p (corresponding to £) for charge Z and mass number A is taken 
as p — A£/Z — 3 GV. For electrons, E in Eq. (16) represents the 
kinetic energy and £ is directly taken as 3 GeV. 

Finally, we take the source spectral index T = 2.2 and the CR 
confinement time within the SNRs as T — 1.2 x 10 5 yrs. Then, 
using the solutions for the CR spectra derived in section 3, we cal- 
culate the CR electron spectrum and the B/C ratio in the Galaxy 
and compare them with the observed data. In Fig. (1), we plot our 
calculated electron spectrum normalized to the data reported by the 
Fermi and the HESS experiments at 50 GeV. This corresponds to 
approximately 0.3% of the total kinetic energy of the supernova 
explosion, taken as Esn = 10 51 ergs, converting into the CR elec- 
trons above 1 GeV. In Fig. (2), we plot the B/C ratio in the Galaxy 
where the solid line represents the ratio expected under the stan- 
dard model which assumes boron production entirely only in the 
ISM, the dotted and the dot-dashed lines represent the results of 
our model which also include additional boron production inside 
the SNRs for rj' = 1 cm~ 3 and 2 cm -3 respectively. In our calcu- 
lation, we consider that boron nuclei are produced due to the spalla- 
tion of carbon and oxygen nuclei and we take the oxygen to carbon 
(O/C) source abundance ratio as 1.4. The data in Fig. (2) includes 
the compilation given in Stephens & Streitmatter 1998 and also the 
recent data from the CREAM experiment (Ahn et al. 2008). From 
the figures, we see that our model explains quite well the observed 
break in the electron spectrum at E ~ 1 TeV without invoking any 
exponential cut-off in the source spectrum and also it explains the 
observed flattening in the B/C ratio above ~ 100 GeV/n. 

Our results can be understood as follows. We see from section 
3. 1 that the spectrum of high energy electrons in the Galaxy follow 
n e (E) oc E^" , "^ 1+fi , where 7 e is the index of electron injection 
spectrum into the Galaxy and f) = (1 — 5)/2 is a small correction 
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HESS 2009 ^ 

Model - 




E (GeV) 

Figure 1. CR electron spectrum at the position of the Earth. The solid line is 
the result of our calculation using Eq. (9) normalized to the data at 50 GeV. 
The data are taken from the Fermi (Abdo et al. 2009) and the HESS (Aha- 
ronian et al. 2008b, 2009) experiments. This normalisation corresponds to 
~ 3 X 10 48 ergs of the supernova explosion energy converting into CR elec- 
trons above 1 GeV. Our model parameters: T = 1.2 X 10 5 yrs, Y = 2.2, 
5 = 0.6, J? = 25 Myr" 1 kpc~ 2 and E S n = 10 51 ergs. Other model pa- 
rameters: wmbr = 0.25 eV cm' 3 , w ph = 0.5 eV cm~ 3 , wpm = 0.2 



eV cm 3 and B 
and in the SNRs. 



6/iG, which are taken to be same both in the Galaxy 



Standard model 




E (GeV/n) 



Figure 2. B/C ratio at the position of the Earth. The solid line is the ratio 
under the standard model where boron is assumed to be produced only in 
the ISM. The results of the present model, where boron production inside 
the sources are included, are shown for two values of matter density inside 
the remnants rf = 1 cm -3 (dotted line) and 2 cm~ 3 (dot-dashed line) re- 



spectively. The ISM density is taken as r\ ■■ 
used are the same as in Fig. (1). 



1 cm 3 . The model parameters 



due to propagation. For the diffusion index 8 = 0.6, we get j3 = 
0.2. In our model, 7 e is the spectrum of CR electrons produced by 
the SNRs (see section 2.1) and is given by 7 e = T for energies 
E < E and 7e = T + 1 for E » E . Therefore, for F = 2.2 we 
get, 



7e 



2.2, 
3.2, 



for E < E 
for E > E Q 



(17) 



which give the corresponding high energy electron spectrum in the 
Galaxy as n e (E) oc E~ A and E~ 4 respectively. The energy at 
which the spectral break occurs depend both on the CR confine- 
ment time within the remnants and also on the total energy density 




E (GeV/n) 

Figure 3. CR path length corresponding to the B/C ratio shown in Fig. 2. 
The dotted line represents the path length in the Galaxy X g and the dot- 
dashed line represents the path length within the SNRs X 3 . The solid line 
represents the total path length (X s + X g ) traversed during their total life 
time in the Galaxy. Model parameters: T = 1.2 X 10 5 , t] = 1 cm -3 , 
r/ = 1 cm~ 3 and S = 0.6. 



Wph of the radiation and the magnetic fields present in the SNRs 
as Eo — l/(aT). For the values of T, w p h and tug adopted in the 
present work, we obtain E = 1.4 TeV. 

The spectra of high energy primary nuclei in the Galaxy fol- 
lows rip(E) oc _E~ 7p ~' 5 , where 7 P is the primary injection spec- 
trum in the Galaxy. For the secondary nuclei, those which are pro- 
duced in the ISM follow a spectrum given by n s (E) oc E~" lv ~ 2& 
in the Galaxy whereas, those produced due to the interaction of the 
primaries within the remnants follow n' s {E) oc E^ Js ^ s , where 7 S 
is the secondary injection spectrum in the Galaxy. We see from sec- 
tion 2.2 that both the primary and the secondary injection spectra in 
the ISM have the same index as the intrinsic source spectrum, i.e. 
7p = 7« = r. Thus, n s is steeper than n' s by E~ s . Therefore, the 
overall s/p ratio follows, 



(n s + n' s ) 



oc — oc E 



for E < E a 



= constant, for E E a 



(18) 



where E a is the energy at which the amount of secondaries pro- 
duced in the ISM equals those produced within the SNRs, i.e. at 
n s — n' s . For our results shown in Fig. (2), we get E a = 7.5 
TeV/n and 2.5 TeV/n for 77' = 1 cm~ 3 and 2 cm~ 3 respectively. 

To better understand the relative production of CR secondaries 
inside the SNRs and in the Galaxy, it is good to compare the amount 
of matter traversed by the CRs during their confinement within the 
remnants to the amount traversed during their propagation in the 
Galaxy. For this, we calculate the path length X as a function of 
s/p as, 



X 



(a p — a s ) 



In 



1 - R s 



(cr s 



(19) 



where m denote the mass of hydrogen atom and R sp represent the 
s/p ratio. Using Eq. (19), we can calculate the CR path length for a 
given R sp . For the path length calculation within the SNRs, R sp is 
given by the ratio N s /N p taken at time t = T (see section 2.2) and 
for the path length in the Galaxy, R sp = n s /n p from section 3.2. 
Their comparison corresponding to the B/C ratio shown in Fig. (2) 
is plotted in Fig. (3) where the dot-dashed line represents the path- 
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length inside the SNRs X s and the dotted line that in the Galaxy 
X g . Also shown in the figure by the solid line is the total path length 
(X s + X g ) traversed by the CRs during their whole lifetime in the 
Galaxy. The average path length within the SNRs is found to be 
X a « 0.09 gm cm -2 for rf — 1 cm~ 3 independent of energy 
whereas the path length in the Galaxy has a maximum value of 
X g « 5 gm cm -2 at around 1 GeV/n which then decreases with 
energy as X g oc E~ s because of the energy dependent escape of 
CRs from the Galaxy. From Fig. 3, we can also see that CRs with 
energies greater than around 7 TeV/n traversed most of the matter 
within the source itself and its effect is seen in the B/C ratio in Fig. 
(2) already at energies around 1 TeV/n. 



5 CONCLUSIONS 

We present a simple model based on the considerations of DSA the- 
ory which can explain both the observed electron spectrum and the 
s/p ratios for a conservative set of model parameters. In our model, 
we assume that CRs after acceleration by SNR shock waves, es- 
cape downstream of the shock and remain confined within the rem- 
nant for some time before they are released into the ISM. During 
this time, CR electrons suffer from radiative energy losses while 
the nuclear species undergo nuclear fragmentations. For a magnetic 
field strength of 6/jG inside the SNRs and assuming a uniform and 
continuous distribution of SNRs in our Galaxy, we find that a CR 
confinement time of 1.2 x 10 yr can produce the observed break in 
the electron spectrum at ~ 1 TeV. Moreover, the hardening in the 
available B/C data above ~ 100 GeV/n can also be explained if we 
assume the averaged matter density inside the SNRs to be ~ 2 cm 
~ 3 . Our results on the B/C ratio which are based on a simple model 
are very similar to those obtained in Berezhko et al. 2003 for the 
case of secondary production inside SNRs calculated for the nor- 
mal ISM density of 1 cm -3 and the CR confinement time of 10 5 yr. 
Their results were calculated using a detailed self consistent model 
of CR production inside SNRs. 

Our model, in its present form, looks similar to the the nested 
leaky box model proposed by Cowsik & Wilson 1973, 1975. How- 
ever, there are major differences in the basic assumptions between 
the two models. They assumed that CRs after acceleration spend 
some time in a cocoon-like region surrounding the sources where 
the primaries interact with the matter and produce secondaries. The 
residence time of CRs inside the cocoon was assumed to be en- 
ergy dependent and after they are released from the source region, 
they undergo an energy independent propagation in the Galaxy. In 
their model, secondary production inside the cocoon dominates at 
lower energies and at higher energies, it is dominated by the pro- 
duction in the ISM. On the other hand, the basic idea of our model 
comes from our understanding of DSA theory inside SNRs. In our 
model, the CR confinement region as well as the confinement time 
are strongly related to the acceleration mechanism itself. More- 
over, secondary production inside the remnant dominates only at 
higher energies while at lower energies, they are dominated mostly 
by those produced in the Galaxy. A more proper treatment of our 
model would be to perform a self consistent calculation of primary 
CR acceleration and their confinement, both of which are strongly 
related to the efficiency of CR scattering around the shocks by mag- 
netic turbulence. Also, we should include the secondaries produced 
during the acceleration process along with those produced during 
the confinement period and if there, also their acceleration by the 
same shock waves which accelerate their primaries. Such a sce- 
nario of secondary acceleration has been considered to explain the 



rise in the positron fraction reported by the PAMELA experiment 
(Blasi 2009). Though it is beyond the scope of this paper to discuss 
this issue, it is worthwhile to mention in relation to our present 
model, that the relative contribution of the accelerated and the non- 
accelerated secondaries to their total spectrum strongly depends on 
the relative time their primaries spend in the acceleration region 
and in the downstream region (see e.g. Kachelriess et al. 2010). 

We emphasize that under our present model, both the high en- 
ergy electron spectrum and the s/p ratio depend strongly on the 
CR confinement time inside the SNRs. Therefore, their data can be 
used to put constraints on the average confinement time provided 
that detailed informations about the magnetic field strengths and the 
matter densities inside the remnants are taken into account. Note 
that the value of the magnetic field strength of and also that of 
the matter density rj' = (1 — 2) cm~ 3 adopted in our present work 
can be different from those expected in SNRs. For SNRs expand- 
ing in an environment consisting of monoatomic ideal gas, their 
values are expected to be equal to the ISM values scaled by some 
constant factor which in the case of strong shocks is approximately 
equal to 4. Moreover, these values can also be different for different 
SNRs and their true values can only be inferred from observations 
of radio or X-ray synchrotron emissions and thermal X-rays from 
SNRs (Bamba et al. 2003, Uchiyama et al. 2007). In future, we will 
include such details in our calculation and also try to extend our 
work to other secondary species like the CR anti-protons and the 
positrons. In addition, we will also include some important aspects 
which we have neglected in our present study like the evolution of 
the remnant, the weakening of the shocks and the energy depen- 
dent escape of CRs from the SNRs. Adding such aspects, particu- 
larly the energy dependent escape, can strongly affect our results. 
For instance, assuming an escape model which follows the same 
energy dependence as in the Galaxy can lead to disappearance of 
the flattening in the s/p ratio at higher energies. 
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